


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1967-12 


Phonon and photon paramagnetic resonance 
lineshapes in ruby 


Crooker, Peter Peirce 


Monterey, California. U.S. Naval Postgraduate School 


http://ndl.handle.net/10945/12825 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
sa Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 
1967 
CROOKER, P. 





UNO MANDO MELEE MUP NT CC IN ay (as 
Seok eames e ta wats Meth 


PRYER PEIRCE CROOKER 







LIBRARY 
NAVAL POSTGRADUATE SCHOGa. 


"10N 














PHONON AND PHOTON PARAMAGNETIC 


RESONANCE LINESHAPES IN RUBY 


by 


Peter Peirce Crooker 
B.S., Oregon State University, 1959 


Submitted in partial fulfillment of the 
requirements for the degree of 


DOCTOR OF PHILOSOPHY 
from the 


NAVAL POSTGRADUATE SCHOOL 
December 1967 


ABSTRACT 


Electron paramagnetic resonance (EPR) and acoustic paramagnetic 
resonance (APR) lineshapes have been measured for four transitions in 
ruby at 5.924 GHz. Cr concentrations varying from .005% to 1% were 
used. Second moments were calculated from the APR lines; half-widths 
were calculated from both EPR and APR lines. A theory is presented 
which yields the second moment in terms of sums -over exchange and 
dipole-dipole interactions between Cr pairs. It is shown that the 
theory predicts moments much larger than measured due to the fact 
that the stronger Cr-Cr interactions lie outside the observed line. 

A qualitative calculation of the moment from the known interactions 
supports this conclusion. Finally, half-widths for both APR and EPR 


lines are compared and shown to be equal. 
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CHAPTER I 


INTRODUCTION 

If a magnetic ion is placed in a crystal lattice such that 
it retains all or part of its original magnetic moment, the 
crystal will display macroscopic magnetic properties. In order 
to understand these properties, it is necessary first to under- 
stand how the magnetic ion interacts with the host lattice, and 
second, how the ions interact with each other. The technique of 
electron paramagnetic resonance (EPR), discovered by Zo TOME Neray 
in 1945, has provided an experimental source from which many of 
the answers can be found. 

In the formative years of EPR, emphasis was placed on the 
location of paramagnetic resonance lines, i.e., the combination 
of microwave frequency, static magnetic field, and crystal 
orientation which would produce resonance absorption of the micro- 
waves. It was obvious, however, that additional information 
could be derived from the lineshape itself, not only with respect 
to its width, but also to its mathematical form, especially as 
it was altered by temperature, crystalline environment, magnetic 
field direction and magnitude, and ion concentration. Impetus was 
given to lineshape study by Van Vleck's classic eee in which, 
rather than attempting the almost impossible problem of calculating 
the complete lineshape, he showed how to obtain the statistical 
moments of the line in a simple and direct fashion and how, knowing 
some of the moments, one could make predictions about the lineshape 
itself. In particular, Van Vleck showed that the line width is due 
to the spin-spin interaction between neighboring magnetic ion pairs 


and that this interaction takes on two forms: 


1) Magnetic dipole-dipole coupling which is caused by one 
ion responding to the magnetic field of another and vice- 
versa. 

2) Exchange coupling resulting from a partial overlap of 
the wave functions of two ions and causing an apparent 
Spin-spin interaction between the ions. 

The result of these mutual perturbations between the ions is to 
take a large number of exactly coinciding energy levels at, say, 


E and spread them out into a collection of levels very near E 


ills T 


Transitions involving E, are consequently spread over a range of 


1 
energies and the line appears broadened. 

Since Van Vleck's paper treated only those ions whose energy 
levels are completely degenerate at zero magnetic field, Pryce and 
Stevalls: were able to extend the treatment to cases where the 
crystalline electric field caused the levels to be initially split. 
In this case, the levels are not equidistant and Van Vleck's approach 
must be modified accordingly. Pryce and Stevens' results are 
especially significant: for their enee, the exchange interaction 
contributes sizeably to the second moment, whereas in the Van Vleck 
case, exchange appears only in the experimentally less accessible 
higher moments. Ishiguro, Kambe, and vem have exploited this fact 
in an analysis of EPR lineshapes in nickel fluosilicate; their 
success has been a strong motivating factor for the present’ study. 
We will develop Pryce and Stevens' result in more detail ina later 
section; suffice us now to say that a measurement of the second 
moment of a line can yield much information concerning the exchange 


forces between neighbors. 


FO 


Late in the 1950's, a new technique appeared for the investi- 
gation of paramagnetic ions in crystals. Called acoustic paramagnetic 
resonance (APR), the method is in many ways similar to EPR. Typically, 
a high frequency (v ~ 10 GHz) acoustic pulse is applied to one end 
of a paramagnetic crystal and allowed to reflect from the other end. 
The height of this return echo is then monitored as a large static 
magnetic field B is slowly varied. At the resonance value of B, the 
echo diminishes sharply; here, the interpretation is that the energy 
levels of the crystal have been shifted by B until there exists a 
pair with energy separation equal to hv for the sound. This particle 
of sound, which is analogous to the photon of EPR, is called a phonon. 
Photons and phonons are not identical in their interaction with ions, 
however, for the energy field-atom interaction is quite different. 
Classically, EPR occurs because the atom has a magnetic dipole which 
responds to the microwave field's magnetic induction H; quantum 
mechanically, this leads to selection rules of Am = 0 or Am = +1, 
depending on the polarizations involved. On the other hand, APR 
is caused by a lattice wave, periodically modulating the crystalline 
electric field environment of the paramagnetic ion. Since the atom's 
wave function is very sensitive to its environment, its orbital 
angular momentum L responds to the lattice and, by the mechanism of 
L-S coupling, the sound wave finds itself ultimately modulating the 
electronic spin S, The perturbation is quadrupolar in nature and, 
as we will show in detail leter, can lead to selection rules of 
Am = 0, +1, +2. Not only is it possible to see different 


— 


transitions with APR than with EPR, but if there is a mixing of the 


m-States, such as when the magnetic field differs from the quanti- 
zation axis, we will see a different blend of transitions between 
the same two energy levels. 

The basic APR technique was developed experimentally by 
Jacobsen, Shiren, and Taine the theoretical basis was developed 
in a series of papers by Altshuler and summarized in a review article 
by Altshuler, Kochelaev, and eetiicae The latter paper also gives 
a summary of the rather limited theoretical status of APR lineshapes. 
For transitions of 4m = + 1, one can show that the second moment is 
exactly equal to that for EPR; but for dm = + 2, the two moments 
are significantly different. For the fourth moments of EPR and APR 
in the Van Vleck case, (i.e., no zero magnetic field splitting), 
Lenton, has shown that they are nearly the same; thus, comparing 
second and fourth moments, the APR line should be wider. No fourth 
moments have been calculated for atoms whose energy levels have 
undergone crystalline field splitting. More work is necessary in 
this direction. 

The fact that APR and EPR lineshapes may differ seems to be 
a little exploited fact. It is of interest to know specifically 
what causes the shape changes and if these changes also occur for 
atoms with crystal field splitting. Meyer, et. cig have reported 
preliminary data on ruby but their study appears to have been 
terminated in the literature. The only other work,to our knowledge, 
is that of ees at the Naval Postgraduate School. His results 


form the initial phase of the present study. 


EZ 


In this paper, we develop a theory for the second moments 
of both acoustic and photon resonance lines, following the 
prescription given by Pryce and stevense We then show how to 
apply the theory to the case in which the magnetic field is at 
an arbitrary angle to the crystal axes and we explicitely evaluate 
second moments for ruby in terms of exchange and magnetic dipole- 
dipole strengths. It will turn out that the second moments for 
EPR and APR lines are identical, but that there is an angular 
dependence resulting from mixing of the energy states when the 
applied magnetic field is not parallel to the crystal c-axis. 

We have also performed experiments on rubies of concentration 
varying from .005% to 1% by weight, measuring both the acoustic 
and photon lineshapes of several transitions. Halfwidths are 
presented for both sets of lines, and moments are calculated from 
the acoustic lines. 

Finally, we compare experimental and theoretical moments, 
discover that the two do not agree, and present an explanation why 
agreement is difficult to expect. We also discuss the halfwidths 


of both sets of lines. 


Bs: 


CHAPTER IT 


GENERAL THEORY 

The technique for obtaining the statistical moments of a 
magnetic resonance line originated with Van eel who showed 
that while the precise shape of the line is prohibitively dif- 
ficult to determine, the invariance of diagonal sums in quantum 
mechanics allows a straightforward calculation of the line's 
moments. The theory was soon expanded in several directions: 
Kittel and Abvahawebe treated the case where the paramagnetic 
ions were dilute, being an impurity added to a host crystal; 
Kambe and cee included the temperature dependence of the 
moments; and Pryce and eer eae derived formulas for ions 
whose energy levels are initially split by a crystalline field. 
More recently, McMillan and fmenowstd have incorporated 
these refinements into one set of formulas. In this section, 
we will rederive the expression for second moments, following 
essentially the approach of Pryce and Stevens. It will be seen 
that the inclusion of phonon induced transitions does not alter 
the theory in any way and that the results apply in a natural 
way to APR. 

We start by supposing a large number of paramagnetic ions 
in a host lattice, such that each ion is characterized by an 
identical set of R energy levels in the absence of any ion-ion 
coupling. Two of these levels, it will be further supposed, are 
separated by energy hv o3 it is the width of the transitions 
between these levels that we wish to calculate. To do this, 


consider the quantum states | of the entire crystal; under 


14 


the assumption of no coupling, these are just products of single 
ion states, so that the energy Ey is characterized mainly by a 
certain number of ions in each single ion energy level. Of course, 
Since there are numerous ways to achieve a given system energy, 

|v) is highly degenerate, containing as it does all possible spin 
combinations yielding the energy Eye Now, if one of the ions in 
|) absorbs a photon or phonon of energy hy os then the crystal 
goes into a new quantum state 18) of energy Eg = E thy | which 

is also highly degenerate. In this model, we therefore see that 
the resonance line is discrete at hy | and no broadening occurs. 

We next introduce a small inter-ion coupling W = 2% W ene 
the summation is over all ion pairs i and j. The effect of W on 
the degenerate states of |v) is to introduce a number of off- 
diagonal terms in the Hamiltonian matrix; we then rediagonalize 
and end up with new states | a) formed from linear combinations of 
the zero order states. The outcome of this process is to lift the 
degeneracy of the zero order states la), so that the energy 
E = EY + {a|Wla). The original energy level has thus been smeared 
out into a very large number of closely spaced levels; our resonance 
line is no longer discrete at hy | but is, instead, an amalgamation 
of transitions from the set of levels | a) near By to the set of 
levels |b) near Bg» all of which appear as a broad transition at 
hy | 

We can now write down the formula for the absorption A(v) 


between states |a) and |b). It is 
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Atv) = (mo. of transitions up) - (no. of transitions down) 


aoe KE. - E_)|<alm|b)|* (hv - E,+ E ) 


b 


-e : (E, - E_)|<alm|b)|* 6[v-(E, - E)/b] 
-E /kT 
we” 


S 


-E /kE -E. /kT 
x . - 


where n. and n, are the populations of levels a and b, \<alm|b) {2 is the 
-E /kT 

matrix element of an arbitrary transition inducing operator M, e is 

the usual Boltzmann factor, the sum in the denominator is the partition 


function over all possible levels E_, and 6[v - (E,- E_)/b] is a delta 


b 


function at hy = EL E.: If it were possible to determine the states 
|a) and |b), we could sum this expression over those states |a) and |b) 
that determine the strength of the line at v; practically speaking, this 
is virtually impossible since then we would have to know all the mixings 
of the zero order states of |) and |B). If we settle for sums over all 
states | a) and |b), however, this amounts to taking the trace of the 
Hamiltonian matrix, which can be done using the zero order states. It is 
this fact that makes moment calculations possible. 


First we make some simplifying approximations with respect to the 


temperature. We write the expression 
io ao 
e -e (E- EB, 


-E,/et -(alW|a) /er -(hv + (b|W[)- ee 
l-e 


as 


e e 


X Chy | + (b|wlb) - (alWla)] 
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Since (b|wW|b) and (alwia) are much less than hv | we drop them 


where they appear as terms and get 
-E /kT - (alWla)/kT - hv _/kT 
Oo O 
e e hy l-e 


A next approximation is to assume (alw|a) << kT and replace 
exp (- (alW |a)/kT) by unity; this is not very restrictive at all 


since (alw |e) is quite small. We are thus left with 


“hy /kT 
hy i & o EL /kI 2 
A(v) = a e | <b Im| a> | 6[v - (E, - E.) pa) 
Ss 
ye 
Ss 


and since the factor in brackets is just a constant, we drop that 
too. In any case, we will ultimately take ratios of the above 
expression so that the bracketed factor will cancel out. 


We now define the lineshape g(v) of the resonance line at ae 


yf kt 2 
a e |<a|m|b) | d[ v - (E, - EL) /h) 


ge(v)dyv = wi 8 
~E/kT Z 
ao [<a|mlb>| 8 v -( EB, - EL) /h J 
a,b 


where the upper summation is taken over just those a,b pairs such 
that hv <(Ey - ES h(v + dv) and the lower summation is over all 
possible pairs. In fact, the denominator is just the area of the 


line. The moments can be defined in terms of g(v) as follows: 
CO 


oe { (v - vo)" g(v)dv 


@) 


2 -E_ /kT Z 
ayia ne : e - E.) |(a|™ |b» | 


area 


=, 
= 


where now the summation is over all pairs. 


= 


Ly 


In order to evaluate this, we first put it in the form of 
a trace; the introduction of projection operators By and 2 baa ile! 


B 
greatly simplify the form. We define P aia 


Py = 3 la) (a| Sum over all | a) in |o) 


so that Py applied to any state selects out only those in the 
manifold |v). These operators will be useful in keeping track of 
the states when we perform the traces. For instance, the area 


can now be written 


~E/kI 
area = XZ e (a|mMlb) (b|M|a) 
a, 
-E /kT 
= eau (a| PMP,M| a) 
a 


and we have a trace. The first moment involves the expression 


-E/ kT P 
X=Z e (E, - E,) |<alm|b)| 
a,b 
~E/ kt 
=Z e ((b|W]b) - (alWla))<al[M|b) <b|m]a) 
aD 
~EV/kI 
7 z e (a | P,MP,WPoM Z P MPGMP Wa), 
and finally the second moment involves 
-E_/kT 
Y= Ze” Ge BE) \alm|p)|* 
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-E_/kT 


@ 
= WP,WP.M + W 
y e (a| PMPs PaWPgM + P WP WP MPGW 


- 2P MP.WP 
Om 


g gMPW la) 


The evaluation of these traces is a tedious but straightforward 
process. We start by invoking one more approximation: that of 
setting exp(-E /kT) equal to unity. Our results are therefore 
restricted to moderately high temperatures, but the calculations 
are considerably simplified. Retention of the exponential factor 
is the basis of the temperature studies mentioned earlier. 

We will use the following notation in what follows: 

M = 2 M. : sum over all ions i 

i 


W=42' Wee : sum over all ions k, £, but not k = &£ 
kk 


| a) = Ja, a.» seme Ja.) = one of R single-ion states 


The expression for |a) is now in the uncoupled representation; 
this is permissable because we are going to calculate a trace. 


First, we calculate the area. 
area = 2% (a|P,MP,M|a) 
a 


Ses Coie ae waar a aie aps M ql@1° “a we ape ae » 
Pqda 


=D 22 ae es a (aa [PM PaM 14.84 a )e°- 
P qa, a, 34 


R 
The sums 2 (a,|a,) are just 2 (1) =R, but the sum over a 


ar, a.=l 
ii iL 


and a is slightly more complicated. Because of Py ja >» must 
q 


be in the Py manifold, so that either ion p or q must be in the 


Ey, 


lower state of our transition. If we call the two levels of the 
transition 1 and 2, then Ja, a Ae ja) » lad or 111); Li@ehis 
were not true, then P, acting back on (aa, | would yield zero. 
Now there are two cases: either p > q or p=q. In the former 


case, we have only the three possibilities 


P P = —_ = 
(ail ae git, | 11) (11[M, | 12) Ghileaye= 0 
P ~ ~ 
(12| P.M, git, | 12) (12|M Pa | 11) 0 
(21|P M aM Pal i ee (21 |, |22) > (21/12) = 0 
so that these terms all vanish. Next we try p = q; we get 


_ 2 
z (ala) i (1| POM Pam | 1) = R <1] | 2)| 


The total area is therefore 


area = 2 Rg! Keay: le a: nr | (1|M |2>|? 
P = / : 


Notice that the form of a is unimportant; all that matters is 
that it connects the two states | 1) and |2). 

We next calculate the second moment, although in less detail 
since the procedure is identical. The expression has three terms, 


the first of which is 


rs 
| 


= % (a|P MP gWPGWPgM| a) 


as aa ti -a,| PM Pal, PW 


a Pe Ke Poo Ja,- -.a) 
i j ki pq a 


Boi 


It turns out that the only non-vanishing possibility here is 


oe ee —_ 
K=% Lf (ala, DB (aa | PM PaW, ,P 


ij a a. a. 
i | | 


‘3 aM tj Pg MJa.a, er 
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where s and t = i or j. Although all the substitutions of i and j 


for s and t must now be considered, we display just one: 


Ns 
Ul 


laps M; PGW, ; POWs j PoM jl) 


(1|M, |2) (21|W, iPoY, | 12) (2|M, ial 


(1|m, | 2) (2|M,| 1) (<21|W, ,|12) (12|W, ,|12) + (21]W; 5[21)(21]w, ,| 12) 


The evaluation of the other possibilities for s and t and of the other 
terms in the moment can be performed in similar fashion. Since there 
are no new assumptions or techniques, we proceed directly to the answer. 
It is 


2 Av") a= = }(11|W, ,}1D) = (21 |W, ,l2) : (21|W, jla2>l° 
| 


2 
+ [¢22|w, ,|22) - (12|W, ,|12) . (2u|w, ,|12)| 


ate 
Maw 


2 
(as |W, ,|4s) oe 2s |W, ,|25)| 
2 
+ [ds |W, .[s1) - (25|W, ,|s2)| 


where we have divided by the area and thereby lost the interaction 
matrix elements. This is the general expression for the second moment. 

Before we close this section, we make a general comment on the 
form of eee If one writes down CE for magnetic dipole interactions, 
one finds a number of diagonal terms and a number of off-diagonal 


terms. The effect of the diagonal terms is to mix the states of Fy 


al 


states among the Pa and other states. In the latter case, we could 
then have transitions between states where E, - Ey = 0, 2hv |; 3hv 
etc., and our moments would include their areas. Since our measure- 
ment is confined to the neighborhood of hv > however, we must incor- 
porate this situation into the theory, namely by truncating a to 
only the diagonal part. This step is not an approximation but 
rather a requirement for obtaining the correct answer. Explicitly 
how the truncation is carried out will be discussed in the next 


section, when we apply our general formula to the specific case of 


ruby. 
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CHAPTER IIT 


APPLICATION TO RUBY 
We next apply our theory to the particular case of ruby. 
We have chosen ruby for several reasons: 

1. The ruby crystal structure and energy spectrum are 
well known. Both energy levels and state coefficients 
have been tabulated as functions of magnetic field angle 
and are readily available in the literature. 

2. Ruby is commercially available in many synthetic forms 
due to its popularity as a laser crystal. It is extremely 
hard and may be polished to a high degree of accuracy-- 
an important point when considering acoustic techniques. 

3. Finally, the ruby lineshape has been a puzzle to many 
workers over the years. We will review this previous 
work in a later section. 


Ruby is aluminum oxide (Al ) in which some of the Al ions 


Z OF 

have been isomorphically replaced by Cr ions in the trivalent state. 

The lattice symmetry is rhombohedral and belongs to the space group 

DS, with four and six sites in the unit cell for aluminum and oxygen 

atoms, respaaaeeeely.. a Figure 3-1 shows the location of the 

aluminum ions in the unit cell; since we will only use the aluminum 

sublattice, we have left out the oxygens for simplicity. The edge 

of the cell is a = 5.128A in length and the apex angle is @ = E33 
A simpler way of looking at this lattice is as a slightly 

distorted hexagonal close packing of oxygens, with aluminums filling 


some of the interstices and the c-axis parallel to the trigonal axis 


of the rhombohedral unit cell. The close packing would be perfect 
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if the rhombohedral apex angle were 53°47' and the trigonal 
coordinates of the aluminum were a/3 instead of .3475 a; neverthe- 
less, we can get a good idea of the lattice if we view it as being 
purely hexagonal. Figure 3-2 shows a gum drop model constructed 

in this spirit, exhibiting the aluminum sublattice both with and 
without the oxygens. The first five neighbors, with which we shall 
concern ourselves later, are also indicated in Fig. 3-2. 

A free Cron ion has a aa electronic configuration and a 4F 3/2 
ground state. The Al site which the Cr ion occupies in ruby is 
surrounded by a slightly distorted octohedron of oxygen ions. Weeks 
has discussed the group theoretical aspects of determining the energy 
levels of the GE ion under the influence of the crystalline field. 
For our purposes, it is sufficient to know that the angular momentum 
of the ion is quenched and an effective spin S = 3/2 describes the 
four resulting levels. These levels are not all degenerate but 
consist of two Kramers doublets: the lower doublet corresponds to 
s. = + 3/2; the upper corresponds to s. = + 1/2 and lies 11.5 GHz 
above the lower ouwie, Upon the application of a magnetic field, 
the doublet degeneracies are resolved and four non-equidistant levels 
appear, but the details of the peal cess are dependent on the magnetic 
field direction. 

The meeting point between theory and experiment in a paramagnetic 
resonance experiment is usually the so-called spin Hemi leone 
we now use this device to. further describe the behavior of Crt in 
Al,0,. Leaving out spin-spin interactions until later, the spin 


Hamiltonian is 


-_ 


= SB i = 2 2 i 
Pig 4 Baiigtis tes SoH) D(S” : S(S + 1)) 
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Fig. 3-2a. Al sublattice of Al,0. (white balls) with surrounding 
oxygen sublattice (gum drops). 





Fig. 3-2b. Al sublattice only. Atoms are approximately in 
horizontal planes. Neighbors of atom O are numbered 1 through 6. 
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where S| and g, are the perpendicular and parallel g-factors, 

B is the Bohr magneton, GH, H,) the applied magnetic field 
components, (S_ oo s) the spin operators, S = 3/2, and D is 
the crystal field parameter. The Vaivecaa of g are pues 1.9840 
and = = 1.9867; the rete of D is 5.75 GHz. Inasmuch as 

g = g we will be abie to use the value g eo 2.00 Pin aie 


alk if 


our calculations with negligible loss in accuracy. Thus, 
=~ = 2 15 
K=g8H .S- DIS) - =) 


where we have substituted 3/2 for S. The axes are rigidly fixed 
in the crystal with the z-axis parallel to the trigonal axis. 

The method of solving this expression for the energy levels 
is to form matrix elements of K with the pure spin states of S~for 
example (m == | 1 | m = 5)-- and from the matrix elements, form 
a 285 + 1 square matrix. If the matrix is already diagonal, such 
as when H ee then the energy levels correspond to the diagonal 
values, one for each pure state. If, however, H and/or Hy is non- 
zero, the Hamiltonian matrix will contain off-diagonal contributions 
from the 5. and EY operators. Then, the matrix must be rediagonalized, 
the pure state vectors will be mixed together, and the energies 
themselves will be different from the previous case. A typical level 
is then described as being in the state | a) = a, |S) t a5|5) a: 


a - >) 7 a, - 3), which is a linear combination of the pure states. 


3 
For the above spin Hamiltonian as applied to ruby, Seenat 
has published values of the energy levels versus magnetic field in 


10° intervals from o to 90°. He has also listed the state coef- 


ficients (the ae in the above example) for each level and the 


2/ 


probabilities for EPR transitions between any pair of levels. 
We have plotted two representative energy level diagrams in 
Fig. 3-3. Figure 3-3a shows the levels when H is parallel to 
the c-axis of the crystal; each level is a pure state as shown. 
Figure 3-3b is for H at 50° to the c-axis; to show the mixing, 
we have explicitely written out one state function. 

Next, we turn to the question of which transitions are 
observable in our resonance experiments. The first condition 
that must be met is that the separation between two levels 
corresponds to hv, where y is the microwave or sound frequency. 
We have plotted the H-9 combinations which will yield resonance 
for our frequency of v = 5.924 GHz in Fig. 3-4. However, even 
if the energy spacing criterion is met, the transition may still 
be unobservable unless we satisfy a second condition: the 
transition strength must be large enough to be observable. 


Transition probabilities are determined from the rule, 


W= - abe, |b)|“e@) 


where ”# = Planck's constant/2m, |a) and |b) are the two levels 
involved in the transition, KH, is the interaction Hamiltonian 
between the perturbing field and the paramagnetic ion, g(E) is 

the lineshape (in units of eee) having area normalized to 
unity, and W is the transition probability per unit time (in 300° Oe 


Lf the interaction has the formiK. =A cos Wt then the transition 


probability becomes 


2 
w= (22) [dalalyl? ew) 
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Fig. 3-4. Values of magnetic field H and magnetic field angle 9 
corresponding to transitions at 5.924 GHz. Transitions used in 
this experiment are shown. 
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Notice that we have also taken g(E) over to g(v); gv) is expected 
to be sharply peaked near the transition frequency vy = (E, - E fh. 
To proceed much further on the subject of transition proba- 
bilities at first seems rather pointless since our study is not 
really concerned with the strengths of the resonance lines. Later, 
however, we will apply our moment theory to the energy levels of 
ruby, and it will be necessary to construct matrix elements almost 
identical to those needed here. Also, if the technique for obtaining 
these matrix elements is developed now, we will be able to show some 
additional features of APR affecting the prosecution of our experiment. 
For this reason we continue, 
Case I, Electron Paramagnetic Resonance (EPR 


For EPR, =-.H.,., where » is the magnetic dipole moment 


int tig 
=—_) 
of the paramagnetic ion and Hy is the magnetic induction of the micro- 
wave field. In our experiment, H, =H, cos Wt in the x-direction 
perpendicular to the ruby c-axis, and since ees gBS then 
coe = - g8H.S_ cos Wt 
int ep Se 

We next form matrix elements witht, .. Since we are only interested 
in relative magnitudes between different transitions, we drop the 


constants and concentrate on expressions of the form |<als by 
xX 


We take two levels 
3 Ll Ll 3 
J) =a) talib tal bral -2 
Ll 
fo) = 512) + bplfy + byl D+ ol = 2 


so that 


S| 


(als |b) = a,b, S18. |5) + a,b 1s |4) + -- 


E 


3 il l 
as = Ts Fy +a, Bea Sale) — 


+ remaining terms. 


The S. matrix elements are found in any quantum mechanics book; in 


this particular case we get 


(als. oy ee a ees tae, + a,b.) + Cab 


al Bef + a,b.) 


pees 
Using the values of a. and b, given by Siegman, this expression 


may be evaluated and squared, thereby giving the strength of the 


transition aie In general, we note that all transitions 


are possible at 8° =e 0° where the states are mixed, but at 9 = Ce 
the selection rule Am ane holds and the transitions 3 = 

nt 3 3 3 : 

= = —~ <5 : 

3 CJ 7? and 9 are forbidden 


Case II, Acoustic Paramagnetic Resonance (APR 

The interaction Hamiltonian for APR has been developed by 
Altshuler, et. a for very general symmetries, and used to 
calculate the sound absorption coefficient 9. We write explicitly 


only the factors which contribute to lineshape: 


6 
c= clE a |(ake |b)|~ sw) 
p=1 | 


where C is a constant, and gone are factors incorporating propagation 
and polarization direction of the acoustic wave relative to the 
crystal axes, and the Ns are collections of spin operators. The 


most general form of this expression has been shown by Domdone to be 


a2 


2 ¥ 
O = a S the ) S 
AL3S— (S + 1)] + B(S,S,+5,5) +B (SS +88) 


Z 


Cs. ce 

A S 
where A, B, C are constants and St = ee en This form has the 
following properties: 


lL) o = 0 tor S = 4, therefore the expression only applies 


to the &>% case. 


2) The A terms have selection rule dm 0 


3) The B terms have selection rule dm = + 1 

4) The C terms produce quadrupole interactions: Am = + 2 
Which of these terms appears in any specific case depends upon the 
polarization and direction of the acoustic wave. For the case of a 
longitudinal sound wave propagating down the c-axis of a trigonal 
crystal, this being our experimental situation, Altshuler's prescription 


yields 


o ="C|(al.Sue fe 1) = 3s°|b)|° g(v) 


where the constants have been absorbed into C. The angular dependence 
of this expression becomes evident when we calculate Oo for our two 


mixed states la) and |b) as before. The result is 


7 7 : 2 
Oo = Cla,b, a a,b) anb, ab, g(v) 


For our purposes, the chief consequence of this calculation is that 

at 640°, all transitions are possible, while at 9 = 0° all transitions 
are forbidden. The latter result is indeed unfortunate, for the 
theoretical expressions for moments would be materially simplified 

if we could perform our experiments on the unmixed states. We do, 


nevertheless, have the proper strengths for off-axis transitions. 
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So far, we have calculated the energy levels of ruby and have 
shown how to derive the transition probabilities of both EPR and APR 
resonances at any angle 9. Our final task is to evaluate, for ruby, 
the second moments of the transitions. 

In the last section, we gave a general formula for the second 
moment of a system of levels, where the number of levels remained 


unspecified. In ruby, there are four levels so that R = 4; hence 


2, lL sy! hs ; 2 
((hdv)") = 7 : }(1a]w, | 11) (21 |W, ,|21) (21|W| 12) | 


2 
+ }(22|W, ,|22) - (12|W, ,| 12) -(21|W, ,|12)| 
2 2 
+](13]w, 113) - (23fm, .]23)|" + [(a3|w, ,|31) - ¢23]W, | 32) 


2 2 
+|<14|w, ,|14) - (24]W, | 24) | + }(14|w, | 41) - (24|W, ,|42)| 


We will use the convention that the levels which are participating 
in the transition are gh and |2), while | 3) and |4) are the other 
levels; since the formula is invariant with respect to interchange 
Om | 1) and | 2) and also with respect to | 3) and | 4), there should 
be no further difficulty in sorting out the proper levels. 


The form of 4 is just that given by Van ae 


— ome IN = Ns 
9 [S..8, - NGGy, Se, 5.) J , a= 
W., = g BO ———_L_+1—_—“1__ +5... 5.55, 
ij 5 1 ia 
Biss 
1j 
: al ‘ : 
where Sa is the distance between Cr ions i and j, and ae is a 
unit vector in Bee E, ,. direction. The first term is,of course, 


the well known interaction between two magnetic dipoles, and the 


last term represents an exchange coupling of magnitude Ay, In 
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the last section, we also discussed the reasons for truncating see 
to include only that part which commutes with the main Hamiltonian. 
This procedure amounts to throwing out all terms of ae which are 
non-diagonal in the ber matrix, and retaining only the diagonal 


terms. Expanding ee in spin variables and performing the truncation, 


we get 
w.. =A_,,f[S>S eee cic +S S_.)] 
ij ij ZZ 2 +i -j -i +] 
Z. 1 
= S Ss -= + 
a oa ie Zi al 6 Sea es 
where 
Ere 
- Le 
ee see Zee 5 Z 
1s 
i 
LJ 


and where Vag is the direction cosine between F, , and the z-axis. 

Since so many quantization conventions are in common usage for 
calculations such as the present one, we should, at this point, 
state explicitly which one we are employing. Our axes are fixed 
with respect to the crystal, with the z-axis in the direction of 
the crystal c-axis. We have chosen the quantization axis to be the 
z-axis, so that when the magnetic field is at some angle to this 
axis, the crystalline and magnetic field directions are "in 
competition" and the levels become linear combinations of the pure 
states (that occur when the magnetic field is along the z-axis). 
The direction cosine Sh occurring in the age is always with respect 
to the quantization direction; therefore, it involves the angle 


between = and the ruby trigonal axis. We thus see that sums such 


as a over all neighbors are independent of magnetic field direction 
i 
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and, being constants of the experiment, can be calculated or 
measured once and for all. 

The next step in obtaining the moments is to form the double 
matrix elements (12|W, ,|12), (12|W, jl2Do-ee sum them, square 
the sums, and add the sums ne il all according to the prescription 
given earlier. First, we put a into the matrix element and obtain 
forms like, say, (12|8,,8_,|12). 
These are easily analyzed; for example: 


(12]8,,8_,]12) = (1B,.]1) (2/s_,]2) 


Ce u 3 ine 3 
= (a, 3] +a Gl + aj(-sl + a,(-$)s, (a,[5) + 4,5) +.a5l-9) + a, |-p 


i” 


ae wil 3 3 i a 3 
X(b, <5 =| +b | ae | +b, Ss (b, 15) +P a>) +b,| >) eye 7) 


*S1s,l5> + aga Gls lot -.. «23s [3 + v,b,Z]s_|a)+ ...) 


iH 


[ 3 (a,a, + a,a,) + 2a,a,][ 3 (bob, + b,b,) + 2b,b,] 


The a's and b's are now looked up in Siegman for the proper magnetic 
field strength and angle, and the matrix element evaluated numerically. 
Double matrix elements can then be formed from the other components 
of Ws in like fashion and the complete moment expression readily 
evaluated. The outcome is of the form 

n>{av") =F, E, [oi + BA, BL, +8 Bei] 
where @, 8, and 6 result from our matrix element calculations. 

In order to calculate the summation over pairs, we must remember 


and incorporate the fact that ruby is a dilute crystal and that only 


a few of the neighboring Al sites of a given Cr ion are occupied by 
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other Cr ions. Kittel and Apeerenicies have shown that E should first run 
over the j occupied sites, giving E = log where k indicates all actual Cr 
neighbours of an occupied site. But this sum is the same as saying that 
the contribution from every site, occupied or unoccupied, is just the con- 


tribution from that site if it were occupied times the occupation proba- 


biIlTty 1, So@enae 


1 


jk k k 
all pairs all occupied all possible 
neighbor sites neighbor sites 
This gives us our final moment 
aa. c gue ~ 2 
as A A 
h’ (Av) 1 Oe eee 


Or 





if we divide A and B by h and express all energies in units of frequency. 


Several conclusion can now be drawn concerning the second moment. 

1. As we noticed in the last section, the transition inducing 
interaction does not appear in our expression and therefore 
the second moment should be the same for APR as for EPR. 

2. The second moment is directly proportional to the concentration 
of the magnetic impurity. 

3. There is an angular dependence in Cissy brought about by the 
dependence of @, 8B, and 6 on the amount of state mixing at 
different angles. If we know the sums in A, and B, (which are 
essentially constants) we can then predict the moment at any angle. 

4. Contrary to the results of roudee we see that for levels with 


crystal field splitting, the exchange does contribute to Cay 


both in 4 squared term and in a cross term with Bis This is a 


a7 


highly favorable result because it removes the necessity for 
measuring the fourth moment in order to get exchange information. 
It turns out that fourth moments, which receive their main con- 
tributions from the wings of the absorption curve where the 
Signal to noise ratio is low, are extremely difficult to measure. 
We now have a choice as regards applying the second moment expression 
to our experimental results. Since a, 8B, and 6 can be calculated from 
matrix elements, and since f is just the known Cr concentration of our 
samples, we could treat the expression as an equation in three unknowns: 


x aS, BAB, and } BL. By measuring three moments with different Qa, 


cod 


B and 6 coefficients, we could then solve the three corresponding equations 


for each sum. On the other hand, since the form of B, is well known, we 


k 
couldvcalculate » pa from a knowledge of the ALO; crystal structure and 
be left with only the remaining two sums, 5 ee and Sf Bie as unknowns. 
We have chosen the latter alternative -- especially since the summation 


Pame Of 2 By has already been performed in a nuclear magnetic resonance 


study by Verber, Mahon, and eamer lia. 
at Se 
on OY, 7 
k = 
k 
so that, dividing by h and squaring, 
2 
2_9 gipt ee a ee 
oae = — ie So remaes |) sista 
k 4 Z k 6 Z 
h ry 


where the summation is to be taken over all possible nearest neighbor 


sites. Verber, et al. quote a value of eae for the geometric factor 


in brackets; using g = 2.00, we thus arrive at the figure of % B as 107 


k 


cHz* ees ence), 
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Z 


2 ~ ne 2 
(Av) =7 [a2 a tin Bnd AB ceeceala,O7 GHz )] 


Without an explicit knowledge of the exchange Be We cannot proceed any 
further theoretically. 

By way of motivation for the present experiment, we should mention 
that the above formula has been used by Ishiguro, Kambe, and ees ina 
study of the second moments of ic <a nickel fluosilicate. By taking 
angles of 0 = 0° and 90°, they showed that B = O and the cross term 


> A,B, does not appear. They then calculated b B ; 


PK r for their crystal 


and solved for the exchange term aod Assuming that a was due 
2 


to the six nearest neighbors only, then EAL = 6K° and they were able 

to pave Because they lacked the cross term, they were not able to 

determine the sign of A, but they were highly successful in determining 
its magnitude. 

If we could achieve acoustic resonance in ruby at Om. if Guris .Oue 
that our cross term would also vanish, and we could solve for EHS 
only. But the transition probability for this angle is, as we have 
shown earlier, zero, thereby forcing us to higher angles where we must 
use the full moment expression. We do have, however, an opportunity to 
determine the sign of the exchange. 

In our experiment, we have chosen four transitions (see Fig. 3-4) 
to measure as a function of Gr concentration: 2-3 Hi at 50°, 2-5 etieeat 
30°, 2-3 Lo at 305. and 1-2 at 80°. Using the state coefficients from 
Siegman, we have calculated aw, B and & by CDC 1604 computer for each 
transition; the results are listed in Table 3~1. APR and EPR lineshapes 


cf these transitions were then taken for six rubies of Cr concentration 


.005%Z, .02%, .05%, .1%, .3%, and 1% by weight. No attempt was made to 
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systematically study angular dependence. 


This completes our application of the theory to the ruby crystal. 
In Chapter VI, when we examine the experimental results, we will return 
to the ideas of this section and attempt to estimate the second moment 
by using the somewhat meager data on A. ,. First, however, we will dis- 


cuss the experiment itself. 


transition 





2-3 Hi, 50° dis 
é O 

2 ene 30 19373 

229 Tice 7.539 

CO 935 


Table 3-1. Coefficients of the second 
moment equation 


re == (a OK 7+ BER B, + 65B,°) 


fomuetour transitions. 
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CHAPTER IV 
THE ACOUSTIC PARAMAGNETIC RESONANCE EXPERIMENT 

In this section we will discuss the acoustic paramagnetic reson- 
ance experiment. In Section 1, we will describe our ruby crystals and 
the bonding of transducers. In Section 2, we discuss the microwave 
equipment, followed by the signal detection scheme in Section 3. 
Finally, we show how to get the second moments in Section 4, and we 
discuss their accuracy. 
1. Crystals and Bonding 

The ruby crystals used in this experiment were purchased from 
Union Carbide Corporation, Linde Division, and Adolph Meller Company. 
All crystals were cylindrical, %"" in diameter by 1 1/8" long, having 
frosted walls, and ends precision ground to a parallelism of 2 arc 
seconds and a flatness of 1/20 wavelength of sodium light. The ori- 
ginal boules were grown by the flame fusion process. Although the rods 
were Specified to have their crystal c-axis parallel to the rod axis, 
we took x-ray patterns and discovered a slight misorientation. Table 


i—(esumna ri Zes elvee informa c 1 On. 


c-axis 
Concentration Misalignment Manufacturer 


Adolph Meller 
Linde 
Linde 
Adolph Meller 
Linde 
Linde 


NO Om fF W BN 





Table 4-1. Data on ruby samples 
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All crystals were of standard quality except for rubies 2 and 5, which 
were Linde Superior Internal Quality for better crystalline uniformity. 

In order to propagate sound pulses through ruby, it is necessary 
to provide a transducing material. We first attempted unsuccessfully 
to bond %"' diameter by %"' long x-cut quartz crystals to the rubies using 
indium as a bonding material as described by aoe. It was, however, 
difficult to control the rather critical bonding temperatures with this 
technique; consequently, our bonds turned out either too thick or pro- 
hibitively nonuniform and we abandoned the method. 

In our next attempt, we tried various organic oils and greases as 
bonding agents, again using our %"' long quartz transducers. The bonds 
were aligned at room temperature in a 500 MHz test apparatus (to be 
described later), but this technique was also unsuccessful. The long 
quartz transducers have two disadvantages: first, they have a low Q and 
hence a low conversion efficiency; secondly, a bond echo is produced when 
the sound wave passes through the bond, thereby wasting energy and often 
obscuring the real echo from the ruby. 

In the final bonding procedure, we used resonant quartz wafers 
instead of the longer quartz transducers. The transducers, purchased 
from Valpey Corporation, were x-cut with a 91 + 2 MHz resonant fre- 
quency and 7/32'' in diameter. The thickness corresponding to this fre- 
quency is about 1.5 mil; in order to guarantee good performance at high 
harmonics, the faces were polished to .00001" parallelism and 1/20 
wavelength (of sodium light) flatness. These crystals are highly re- 
sonant, even at 6 GHz,and the conversion efficiency is much higher than 


with the larger transducers. Also, there is no bond echo. 
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We give a step-by-step description of the crystal bonding process 
in the Appendix; suffice it to say here that the bonding agent is Nonag 
Stopcock Grease (Fisher Scientific Co.) and that the transducers can be 
successfully fixed to the end of the ruby. Fig. 4-1 shows a completed 
bond. 

After making the bond, it was checked out on a 500 MHz pulser- 
receiver setup to see if it was good. The checkout apparatus was very 
simple: a 1 watt General Radio 1209 BL unit oscillator, receiving 
trigger pulses from a Teletronix PG-200AA pulse generator, fed 500 MHz, 
1M@sec pulses directly to a Navy AN/APR-4Y receiver. An adjustable 
frequency reentrant cys and a tuning stub were then coupled in by 
tee fittings. After inserting the quartz-ruby bond into the test cavity, 
the resulting receiver echoes could be directly observed on an oscillo- 
scope. Over 100 observable echoes resulted, depending upon the bond 
quality; therefore, poor bonds could be quickly detected and redone. 
Typicable bonds, as could be seen by observing interference fringes, 
had a thickness of a few wavelengths of light, but the bond thickness 
was not particularly uniform. This fact did not seem to affect the 
high frequency echo-train too much; indeed, the lack of rigidity of the 
thin transducer may have helped the bond to remain intact at low tem- 
peratures. 

2. The Microwave Equipment 

A diagram of the J band hypersonic spectrometer is shown in 
Fig. 4-2. A Teletronix PG=200AA@pulise generator prop@ded [igsec, 

100 V pulses at a repetition rate of 103 Hz to a modulator from a Navy 
SU radar. The modulator then drove a watercooled Litton L-3467 magne- 


tron, which in turn produced 5.924 GHz pulses of 1.8 kw peak power. An 
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attenuator was used to reduce the power during the experiment to 100 
watt or less, depending upon the quality of the bond. In order to re- 
duce heat leakage into the helium bath, the signal was coupled from the 
waveguide into a copper plated, thin wall, stainless steel, coaxial 
waveguide, to which the sample cavity was attached. 

The cavity was of the reentrant type (Fig. 4-1) with internal 
dimensions of 1.63 cm diameter by .712 cm depth. A raised post of 
height 1.25 mm served to direct the electric field into the transducer. 
The microwave power was coupled into the cavity by a small loop oriented 
with its plane perpendicular to the cavity magnetic field. In order to 
adjust the resonant frequency of the cavity to a transducer harmonic, 
the ruby was mounted in a collar; the position of the collar and, hence, 
the depth of the ruby in the cavity could be controlled while the cavity 
was immersed in liquid helium. Prior to the experiment, but after cool- 
down, the ruby position was easily adjusted for maximum echo height. 

The cavity Q was about 400. 

The received echoes were passed through a TR (transmit-receive) tube, 
which provided a dead short to the main pulse and protected the receiver, 
to an LEL CAC-2 mixer-preamplifier. A Varian X-26E klystron provided a 
beat signal 30 MHz away from the main frequency. The 30 MHz difference 
signal was then amplified by an LEL 3378 IF amplifier which put out the 
rf signal and the detected envelope. The 30 MHz rf signal was used 
for local oscillator control; it was fed to an LEL 3379 discriminator 
which then provided a correction voltage for the reference klystron. 

The detected signal from the IF amplifier was then available for oscil- 
loscope observation or, as we shall describe in the next section, for 


signal processing. Figure 4-3 shows a typical echo train, first when 
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(a) 


(b) 
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the magnetic field is off resonance, and then when it is on resonance. 
3. Pulse Detection and Calibration 

In order to monitor the resonance, the echo train was passed 
through a gate which selected one pulse for observation. Since the 
duty cycle of this pulse was so low (107*) it was necessary to stretch 
the pulse out; this was done very simply by passing it through a diode 
to a large grounded capacitor. The capacitor had a rapid charging time 
(< 1 sec) due to the low forward resistance of the diode and a long 
discharge time determined by the backward resistance of the diode and/ 
or the input impedance of the instrument detecting the capacitor's 
charge. In this way, the pulse was converted from a 1 Msec spike to an 
exponential decay with decay time T™ 5 msec. This signal was then 
added to a positive going 103 Hz square wave such that the signal 
appeared during the first half of the repetition period and the square 
wave appeared during the second half (when the exponential tail of the 
signal was small). The square wave thus served as a comparison for the 
signal, converting it to an AC signal whose average value off resonance 
was zero. Detection of this signal was then effected by an EMC model RJB 
phase sensitive amplifier, the output of which provided a DC voltage for 
a Varian G-10 strip chart recorder. Synchronization signals for the 
modulator pulser, gate pulser, square wave generator and phase sensitive 
detector were supplied by a Hewlett Packard 200 AB sine wave generator. 

The magnetic field was provided by a Varian V-4012-3B 12" magnet 
operated from a V-2100A power supply and controlled by a Varian V-FR2100 
Fieldial control unit. Uniformity of the field was better than .05 G 
over the sample, and the Fieldial control unit gave extremely good 


stability in time. Calibration of the field was obtained by using an 
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Alpha AL67 NMR gaussmeter, modified so that the NMR frequency could be 
extracted, amplified, and read on a Beckman 6127 frequency counter. 
Since the NMR frequency was stable to™ .Ol KHz, the accuracy was about 
one part in 10° for a 10 MHz NMR vyoltaze. 

Since quartz transducers are highly absorptive to high frequency 
sound at high temperature, the experiment was performed in liquid heliun. 
The low temperature apparatus was of the standard type: an outer dewar 
filled with liquid nitrogen surrounded an inner dewar filled with liquid 
helium. The vacuum jacket on the inner dewar was connected to a small 
forepump only; when liquid helium was introduced this jacket was sealed 
off and allowed to cryopump. Control of the inner jacket vacuum turns 
out to be essential in crder to achieve slow cooldown times and prevent 
breaking of the transducer bond. 

In order to reduce bubbling in the cavity, the helium was pumped 
below the A-point. The temperature at which all lines were measured 
was 1.7 K. 

A typical APR lineshape was calibrated as follows. The magneto- 
meter was set to some field value, then, when the magnetic field swept 
through the magnetometer resonance, a pen mark was put on the chart 
paper and the frequency was read. The magnetometer was then reset and 
the process repeated, providing for about six cor seven marks per chart. 
The frequencies were then converted to magnetic field values and plotted 
against chart distance; in all cases, the relationship was linear. The 
accuracy of this procedure was limited maily by the rate at which the 
magnetometer signal crossed the oscilloscope face; accuracy is gener- 


ally + .5 gauss and in no case worse than + 1 gauss. 
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Immediately upon completion of a chart trace, a signal calibration 
was made by setting in known amounts of attenuation with the precision 
attenuator and recording the corresponding recorder pen positions. This 
procedure accounted for all non-linearities introduced between the 
attenuator and the recorder--virtually the whole detection system. A 
calibration curve of dB attenuation versus pen position was then made 
and the chart ordinates corrected to read dB. The accuracy here depended 
upon the amount of background noise, but generally can be set at + 5%. 
It was necessary to avoid sweeping the resonance too fast on some reson- 
ances, since the trace can be controlled by the lock-in amplifier inte- 
gration time rather than the signal itself. This effect was calculated 
beforehand and a set of maximum allowable chart slopes was obtained and 
never exceeded in the experiment. 

The overall sensitivity was a few tenths of a decible, however 
on the lowest concentration ruby (.005%) no signal was detectable. 
Figure 4-4 shows two sample curves, one very good and the other rather 
poor. These graphs represent the extremes of our data. 

4. Data Evaluation 

Once the data was taken, the next step was to obtain the second 
moments. Each chart was first digitized to between 40 and 80 points 
and then analyzed by IBM 360 computer. The computer program performed 
the following calculations: 

l. For a given transition, five values of the energy level sepa- 
ration v corresponding to various magnetic fields H were 
entered. For example, in the transition 2-3 Hi at 50° which 
occurs at H = 3000G, v was entered for H = 2000G, 2500G, 


3000G, 3500G, and 4000G. The program then fit these pairs 
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with a five degree polynomial and generated a frequency point 
for each digitized value of H. This process converts the 
magnetic fields to frequencies, automatically accounting for 
asymmetry caused by nonlinearities in the v-H curve. 

2. The area of the line was calculated using trapezoidal incre- 
ments. Each ordinate was then divided by the area, thereby 
normalizing the line and yielding the true lineshape g(v). 

3. The second moment was calculated. 

4, Values of H, v, and g(v) were read out along with the value 
of the second moment. 

It must be remembered that the moment calculation is highly sensi- 
tive to the behavior of the line's wings, yet it is in the wings where 
themsignal to noise eatioem@eue poorest. The limiting@actonment our 
experimental accuracy is just this effect, rather than the approximations 
inherent in digital integration. We therefore feel that our moments have 
an accuracy ranging from + 10% for the strong lines, up to + 504 for 
weak, and consequently noisy, lines. 

In addition to the moments, the widths of the lines at the half- 
power points were obtained by interpolating between the computer print- 
out values. These half-widths were dependent mainly on the frequency 


calibration and therefore are felt to be accurate to a few percent. 
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CHAPTER V 
THE ELECTRON PARAMAGNETIC RESONANCE EXPERIMENT 

Of the equipment described for the APR experiment, only the 
crystals and DC magnet system are the same for EPR. The microwave 
equipment and detection system for EPR are different, as we will 
describe in the first part of this chapter. The second part will 
be devoted to the data analysis and experimental accuracy. 

1. Microwave Equipment and Detection 

The J band microwave spectrometer used in this experiment is of 
standard superheterodyne design as shown in Fig. 5-l. Both ae and 
local oscillators were Varian X-26E klystrons operated in the CW mode 
at 5.924 GHz. The signal was first fed through an E-H tuner, isolator 
and two 20 dB attenuators to a hybrid tee. Although the hybrid tee was 
initially set up to work as a microwave bridge, this proved to be un- 
necessary in the experiment, with the result that 50 dB of attenuation 
was permanently introduced into the reference arm. 

The signal side of the hybrid tee was connected first to a fre- 
quency meter, then to a directional coupler, and finally to the cavity. 
Having the frequency meter in the cavity arm allowed us to lock on to 
the meter resonance if necessary, but this option was never used. The 
directional coupler was useful for tuning the klystron to the cavity and 
measuring the frequency, neither of which measurements required ,the re- 
Solvers 

The cavity was a section of J-band waveguide ‘7.18 cm long operating 
in the TE mode. The ruby was inserted through the shortest dimensions 


G2 


of the cavity midway between the aperture and the cavity end; the trans- 
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verse magnetic field was therefore perpendicular to the ruby c-axis. A 
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small metal or plastic screw could be inserted in ie cavity bottom for 
fine tuning; coarse tuning was effected by changing the aperture size, 
normally set near, but not on, critical coupling. Cavity Q's were about 
3000 with the ruby in place. 

The return signal from the cavity was sent directly to an LEL CAC-2 
mixer-preamplifier and the derived 30 MHz difference signal amplified 
in an LEL 3378 IF amplifier. The rf output from the amplifier was used 
two ways: first, it was monitored on an oscilloscope so that its change 
during resonance absorption could be checked; secondly, it was fed to an 
LEL 32B93 discriminator (modified for CW operation), the output of which 
was used as a correction voltage to lock the reference oscillator to the 
signal. A detected rf envelope, also provided by the IF amplifier, was 
fed directly to an EMC model RJB lock-in amplifier, the DC output of 
which then drove a Honeywell Electronik 19 recorder. 

One feature of the spectrometer is not shown in Fig. 5-1. In order 
to protect against frequency drift and, more critically in a lineshape 
measurement, to guarantee that only the absorptive part of the complex 
susceptibility was observed, the signal klystron was locked onto the 
cavity resonance. To do this, the signal klystron reflector voltage 
was modulated at 10 kHz by an external sine wave generator. Due to the 
cavity characteristic, the received signal exhibited this modulation 
as an amplitude modulation, with the following properties: If the 
klystron frequency was situated on one side of the cavity, an approxi- 
mate sine wave resulted; if the frequency lay on the other side, a sine 
wave with the opposite phase appeared; if the frequency coincided with 
the cavity center frequency, only a small second harmonic appeared. The 


) 


Sine wave was then phase sensitively detected and the resulting error 
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voltage added to the signal klystron reflector such that the signal fre- 
quency followed the cavity while the cavity shifted during passage 
through the line. A considerable reduction in the signal to noise ratio 
also resulted, indicating a lack of frequency stability without the 
automatic frequency control (AFC) feature. 

The static magnetic field was provided by the same system used in 
the APR experiment. This field was modulated, however, by passing a 
sinusoidal current through coils wrapped about the magnet pole faces. 
The modulation frequency was 103 Hz, a frequency which is easily 
separable from the 10 kHz AFC signal but slow enough to accommodate the 
long relaxation time of ruby. The effect of the modulation was to cause 
the cavity reflection to oscillate correspondingly as determined by the 
shape of the ruby absorption characteristic. Since it was only the 
103 Hz component of the cavity reflection that was ultimately seen by 
the lock-in amplifier, the resulting chart was the derivative of the 
true absorption line. In order to avoid distortion of the line, the 
amplitude of the modulating field was never more than 3G peak-to-peak. 

The low temperature apparatus was not used for the EPR experiment 
since the lines could be easily observed at room temperature. 

Most of the curves were obtained with about 10 LW of power to the 
cavity, although this value was raised to as much as 100 LW in some 
instances. Several curves exhibited consistent anomalies in their 
lineshapes, perhaps due to some paramagnetic impurity in the crystal. 
These were largely unexplainable; in any case, data from these curves 
was omitted in the final results. A typical derivative trace is shown 
in Fig. 3-2. 


We now discuss several points relative to obtaining good EPR data. 
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Sample EPR derivative trace. 
magnetic field calibration. 
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In order to observe the true lineshape of an EPR transition, it is 
necessary to measure only the power absorption of the sample, but this 
quantity is not, in general, simply related to the cavity reflection. 
it En and EB sue are the signals into and away from the cavity, R is the 


on-resonance effective cavity resistance, and Z is the characteristic 
O 


impedance of the waveguide, then 





E R- 2 
out _ O 

Eo R. sez 
Ln O 


and we see that Eats does not vary linearly with R. The resolution of 
this problem has been discussed thoroughly by ener ae he has shown 
that in the limit of small (ideally, differential) changes in the cavity 


, AE is proportional to the 


durine@enesonance. ste wchancewin 
Q g g ; ae 


We 
true lineshape. We have assessed the errors involved in this technique 


with the following result: if E changes by 1% during the resonance, 


out 
then the error in the true lineshape is also 1%. For this reason, the 
full rf envelope from the receiver was observed during every run. If 
the height of the envelope changed by more than » 1%, the cavity Q was 
lowered by changing the cavity iris size, and this lessened the change 
mE due to the ruby resonance. 

out 

Another difficulty that occurs in EPR is the fact that the data 
appears in derivative form. In order to obtain the moments, one must 
first integrate the derivative to obtain the absorption curve, then 
integrate again in order to finally obtain the area and moments. The 
effect of the first integration on the moments can be very critical, for 


if the two lobes of the derivative curve differ in area by only a few 


percent, the absorption curve will not return exactly to the baseline 
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in one of its wings, and it is mostly the shape of the wings that de- 
termines the size of the moment. 
In our experiment, it turned out to be very difficult to match 
the lobe sizes of the derivative curves to better than a few percent; 
for this reason, it was not possible to get moments. We feel that with 
more time and considerable effort, the accuracy of the raw data could 
have been adequately improved; however, in light of the fact that the 
APR moments resisted quantitative interpretation anyway, we restricted 
our results to EPR linewidths, for which the present data was sufficient. 
2. Data Analysis 
The analysis of the EPR data was identical to the APR analysis 
except in the initial conversion of the line to an absorption line. 
Values of H were derived from gaussmeter frequency markers on the chart 
and plotted against chart inches. The resulting calibration curve was 
linear in every case, yielding magnetic field values for each point on 
the chart trace. The chart was then digitized to between 50 and 100 
points and analyzed by IBM 360 computer. The computer program performed 
the following functions: 
1. The derivative trace was integrated using trapezoidal incre- 
ments and magnetic field abscissas. 
2. The magnetic field was converted to frequency values v, in 
exactly the same way as in the APR analysis. 
3. The area of the absorption curve was calculated and used to 
normalize the absorption line ordinates so that the area was 
unity. The ordinates thus became values of g(v). 
4. Moments were calculated. 
5. The program printed out g(v) versus v and H for each digitized 


pormt, 


ay 


Half-widths, corresponding to the full width at half power, were 
obtained by interpolating between the computer print out values. 


The half-widths are thus felt to be accurate to within a few percent. 


60 


GHAPTER VI 
RESULTS AND DISCUSSION 

in this section, we will present the measured second moments for 
APR and relate them to theory. Near the end of the section, we will 
also display EPR and APR half-widths and discuss them qualitatively. 

The second moments of our four transitions are plotted in Fig. 6-1 
as a function of Cr concentration. We note several features of these 
curves: 

l. At 1% concentration, the curves have widely different moments, 

varying from .007 aie for. 7 3 GHz? 

2. The curves approach a constant value of the order of tom He 
as the concentration is lowered below about .05%. The effect 
1s independent of the transition. 

3. From the standpoint of our theory, the most important result 
is the extreme smallness of the measured moments. Using only 
the dipolar interaction, the second moment is predicted to be 
(ey = (concentration) (EB) ~ (17) CLO? GHeay 1 CHz* for 
a 1% Cr concentration ruby. But as experimentally measured, 
the largest moment at 1% is ~ .07 cH” 

We begin our discussion of the data with an explanation of the moment 
size. 

The fact that measured ruby second moments might be smaller than 
theory predicts has been alluded to by Manenkov and Reaetoy and 
Grant and Sicrandbereaee although neither of these authors actually 
measured the moments. Manenkov and Fedorov measured EPR line widths 
€ at © = O° and assumed either a Gaussian line, where RO = sien 


or a Lorentzian line with cutoff at 5e (a Lorentzian lineshape must be 
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terminated somewhere in the wings or Ges diverges), where 

: 2 
(Av ) = 1.35€°. They then showed that neither of these moments is 
large enough to account for the dipolar part of the theoretical value. 
In point of fact, however, the line is neither Gaussian nor LorentZian, 
so we can consider their values of (ay) to be just an estimate. 

Grant and Strandberg have done a study on line shape in ruby and 
have inferred most of the second moment is not present in the observable 
line. Our measurements indicate that this is indeed the case. [In the 
following, we will show how to estimate the size of the interactions 
involved in a resonance line; then, we will compare the measured widths 
with theory; and finally, we will try to calculate the second moment by 
summing only the interactions from neighbors outside a certain radius. 

We first demonstrate that the overall width of a resonance line, 
from where the line is first measurable to where the line returns to zero, 
is of the order of the largest ion-ion interactions that are causing the 
shape. For example, consider two identical uncoupled spin % particles 
in a coupled representation. There will be three energy levels: one at 
+E corresponding to both spins up (++), one at -E corresponding to both 
spins down (--), and two degenerate levels at zero energy corresponding 
to up-down (+-) and down-up (-+). Only one transition of energy E is 
observed from this system, if we discount double spin flips. 

We next ask the effect of turning on a small amount of coupling 
between these particles, say of the form SOI 59° Since the Hamiltonian 
matrix is still diagonal, the change in energy can be found simply by 


forming matrix elements of the form (a|s_ Sola) and noticing that while 


il 
the (++) and (--) states both move up a small amount A, the degenerate 


states move down by A. There are now two transitions of size E + 2A and 
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E - 24, respectively; our original transition has "broadened" into a 
line of overall width 4A. 

One can also suppose a small off-diagonal coupling of the form 
S41 no te 5 S45 between our original uncoupled particles. In this 
case, we must rediagonalize our Hamiltonian since the perturbation is 
between degenerate levels, but the end result is similar: the degene- 
rate levels are mixed, one moving up by A, the other moving down by A. 
The original transition at E now consists of two lines at E + A and 
E - 4+ they width is thuseZ4, 

Of course, this simple analysis cannot be taken too far. It is 
not the interaction of just ions i and j that causes broadening ina 
large crystal; it is one of the states composed of all the mutually 
interacting ions of the crystal that yields a resonant frequency dif- 
ferent from the single ion frequency. Consequently, when we say that 
the interaction between ions i and j appears Ay away from the main 
resonance frequency, We can expect only order of magnitude accuracy. 

We define, then, the "full-width" of a line to be the distance 
in frequency from where the signal is first detectable above the noise 
to the center of the line. We also assume that if two Cr ions interact 
with an energy €, they will produce a displacedline at frequency 
Avy = ¢/h away from the single ion frequency. In this sense, a measure 
of the full-width of an observed line gives us a criterion for deter- 
mining which interactions are contributing to the observable line. 

The full-widths of our lines are given in Table 6-1; we note that 
the full-widths vary from 1135 MHz down to 85 MHz. In order to compare 
these values to actual interactions, we calculate first the dipolar, 


and then the exchange, contributions to the second moment. 
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Table 6-2 is a summary of the dipolar contribution BL from each 


nearest neighbor for the first 20 nearest neighbors. It is obvious 

that the individual contributions of many of the neighbors is more than 

the total observed second moment itself. 

ik has been given by Statz, 
22 4 , 

et al. ; they quote a value of (1.17 + .15) x 10 GHz for the first 


An estimate of the size of the exchange A 


neighbor at 2.73A and about 15 GHz at 5.73A (eleventh nearest neighbor). 


The contributions fo eA. are then 1.37 x 10° capa and 225 ai 


k 
respectively, for these ions. The contrast between theory and experi- 
ment here is even more drastic than with the dipolar case, the differ- 
ence being many orders of magnitude. 


Ruby Number 


Z 5 4 5 6 





223 Hi, 50m 620 
2-3 Hie 308 1135 
Qed Lonmeen = 
12a Oe 880 


Table 6-1. Full-widths of APR lines. 
Widths are from beginning 
of line to center. 
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Neighbor Number r, A 9, deg B-, CHS 


Ke? 
0 it -- -- =~ 

L i 2a 8 180 Sicko 

Z 5) 2.804 LO... 6 9.48 

3 3 3,180 Be aL Mee, 
4 3 3.498 3) te) s Oe 
4! 3 3.498 LZ2GR2 nO72 
3 1 5.67, 0 8.49 
6 3 4.758 90 goul® 
Go" 3 4.758 90 LO 


Table 6-2. Dipolar contribution of the second 
moment from each of nearest 20 
neighbors. (After ref. 14 and 24). 
The sum of By over all neighbors is 
107 GHz. 

The picture that emerges is now fairly clear. The theoretical 
second moment as we have calculated it is correct, but the contribution 
of the near neighbors to this moment is so immense that the transitions 
caused by these interactions lie in the far wings of the line and are 
unmeasured. For instance, the nearest neighbor, which lies only 2.73A 
away from the origin on the z-axis, has an exchange contribution which 
ii 
is 10 times bigger than the observed moment; furthermore, its dipolar 
contribution is, in itself, 10 times larger than the measured value. 
Indeed, for the first few nearest neighbors, the exchange coupling is 
so strong that the ions form pairs and tend to produce pair spectra 
which are quite independent of the single ion line. There is, in fact, 
a considerable literature built up by workers trying to locate the 

25-27 


spectra of these pairs, not without some success . For our purpose, 


however, which is to explain the observed second moments, it is obvious 
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that we must exclude such pairs from the summations over neighbors. 

In order to get a quantitative understanding of the observed line, 
we are faced with a dilemma. We know we must exclude the near neighbors 
in our summations and sum only outside some radius o> also, to avoid 
circular reasoning, we must choose tO without reference to the second 
moment. If, instead of calculating the second moment (which is insensi- 
tive to line structure), we were able to calculate the exact line shape 
by knowing which atoms contribute to each part of the line, our approach 
would be obvious. Lacking this straight-forward situation, we resort 
to the following technique. 

We assume that if a Cr ion interacts with a neighbor at radius 
r such that the interaction strength is greater than the full-width 
of the line, then that interaction is not being observed and must, 
consequently, be excluded from our summation. By matching the full- 
width of a line to the interaction size, we can therefore arrive at a 
cutoff radius for each term in our second moment. 

Obviously, this argument is not without a certain amount of cir- 
cularity; we use it only as a last resort for determining te We can 
therefore expect no better than order of magnitude agreement between 
our measured and calculated value, especially insofar as we will have 
to make certain assumptions about the relatively unknown exchange 
energies. 

We will choose two radii: a is that radius outside of which we 
will sum the exchange contribution, and r, is that radius outside of 


d 


which we will sum the dipole contribution. Since we and ty will turn 
out to be fairly large, we will perform the sums by integration outside 


of these radii. We start with the exchange part. 


A7 


Since the data on exchange is meager, we have fitted the pre- 

viously quoted two values of Statz, et See the equation 
toeaalnarse 

where C = 9.93 x 107 CHz-A?, and n = 9.00. This is, we should point 
out, a very arbitrary choice for A at all radii. Mollenauer and 
saice have measured A for the first four nearest neighbors; their 
values do not obey an r" law. But even a knowledge of A(r) for near 
neighbors does Little to help us predict A for the far neighbors. The 
radial dependence of A, for instance, depends critically on how many 
oxygen ions lie between the Cr ion pair; also, one can obtain both 
positive or negative values for A, indicating ferromagnetic or anti- 
ferromagnetic coupling, respectively. (Ke, however, is always positive.) 

We now choose the transition for our calculation: 2-3 Hi at 50° in 
sample 6 (1% Cr concentration). This line has a full-width of 620 MHz 
(Table 6-1) and a second moment (iy = ,036 GHz? (Fig. 6-1). The ex- 
change cutoff radius t > obtained by setting [Al = 620 MHz and solving 
for r, is 8.1l6A. This is a large radius-- asphere of this radius con- 
tains some 40 neighbors-- so we will approximate the summation outside 


a by the integral 


e 2 ae 2 
eo § A® 9 4nr“dr 
r r 
e e 
where § = number of ions per unit volume = .0574 Aone and A = C/r™ 


as before. The result is 10.1 Cz". 


The next sum is the dipole sum. We recall from Chapter II that 


the dipole-dipole energy for neighbor k is 


-39°B2(3y* - 1)/2 


: 5) 
k 


By = 
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Ignoring the directional factor and setting B, = 620 MHz, we can solve 


k 
for the dipole cutoff radius ty and get a 6.30A. The summation 
9 | 

over B. then reintroduces the directional factor: 

2m TT oo 

2 BL = i ao \ aime 20 | 0 BL aa 
O O ty 
2 eee : J 

Theuwfactormi@y2 ear +) is thus averaged over all directions and 


we get b BL - 4. 63GHz~. 


The last Ssumeus ye 


Hee which vanishes because the angular 


average of 3/2 Ve - 4 is zero. 


The moment expression is, then, 
a oe ee 2 
(Av™) = 7 (@ ZA, + 6 ZB") 


Inserting f = .01, @ = 9.47, and 6 = 1.37 for this transition (see 
Table 2-1) we arrive at (iy) = ,255 Che as an estimate of the second 
moment. Although our measured moment of .036 CHa is some / times 
smaller than this, we feel that we have arrived at order of magnitude 
agreement with experiment. 

Another feature of the second moments of high concentration is 
their variability trom transition to transitiong Sin view oreear 
inability to accurately predict these moments, it is difficult to 
assess the cause of this variability. If we assume that DAL 3 


py AB and »& By are constant for all lines, as theory predicts, then 


we should be able to solve for the sums using three moment equations 

and the knowcoefficients. Hopefully, the angular dependence would 

then be due to the angular variation of @, B, and 6. We have attempted 

this fit unsuccessfully; in fact, we have even obtained large negative, 
Z 


and therefore uninterpretable, values for & Be : 
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We could postulate that variation is caused by a wandering of 
the crystal c-axis as we proceed along the sample rod axis. Wenzel 
and aioe have discussed this effect in order to explain the low- 
concentration linewidth; in our case, it will affect those transi- 
tions whose angular variation dH/d® is large. From Fig. 2-4, we 
see that the transitions 1-2 at 80° and 2-3 Hi at 50° are most eligble 
for this type of broadening. We have calculated that, for a half- 
width of 300 MHz, the c-axis would only have to wander a degree or so, 
which for a high concentration ruby (1%) is not at all unreasonable. 
We feel, however, that while these transitions are indeed the larger 
observed ones at 1%, a careful angular study would have to be made to 
confirm this explanation. 

The third feature of the second moment curves is the residual 
moment of .002 GHz? at low concentration. This effect has been ob- 
served as a residual half-width in EPR measurements, and Grant and 
Seerisere have reviewed the various explanations. More recently, 
Wenzel and anes have had some success on the problem by using various 
models for crystalline imperfections. A calculation of the EPR second 
moments at low concentration has been made by Laurance, McIrvine and 
manner who assume that the residual width is due to hyperfine inter- 
action with the neighboring aluminum nuclei. While our second moments 
are for APR, their value of .000029 GHz" is orders of magnitude too 
small. 

Inasmuch as the experimental second moments for APR have resisted 
quantitative interpretation, we confine our EPR results to the half- 
widths of the transitions only. (By half-width is meant the full 


width of the line at the half power point.) Figure 6-2 gives the EPR 
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half-width; Fig. 6-3 presents the corresponding APR half-widths. The 
curves are quite similar. Both APR and EPR half-widths are about 

320 MHz at 1% concentration and they both decrease to approximately 

70 MHz at .05% concentration and lower. The angular dependence of 

the lines is qualitatively the same: almost no angular dependence at 
1%, while at low concentration the ordering in magnitude of the tran- 
sitions is the same. Table 6-3 compares the half-widths at high and 
low concentration. We feel that ane differences that occur are quite 
minor, but that a systematic study of linewidths versus angle would be 
useful to confirm this in detail. As a further complication, it should 
also be borne in mind that the half-width measures the central part of 
the line, so that while the second moment deals with the distant 
neighbors of a Cr ion, the half-width is caused by neighbors even 
further away. A theory of half-widths is therefore even more dif- 


FiCUuLt to ConsErdce. 


ee ae 
Transition EPR, MHz APR, MHz EPR, MHz APR, MHz 


2-3 Hi, 50° 308 
2-3 Hi, 30° 56 56 308 315 
2-3 to 30 | 7B 91 470 sae 
1-2 at 80° 72 67 S71 336 


Table 6-3. Half-widths of EPR and APR lines 
at two concentrations. Half-widths 
are the full widths at half-power 
points. 
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We now finish this ection with a summary: 


ie 


We have rederived the general second moment formula of Pryce 
and Stevens and have shown that it is also applicable to APR. 
We have then applied this theory to ruby and arrived at an 
expression for (iy 7) in terms of angular dependent matrix 
elements, Cr concentration, and sums over A, and Bie the 
exchange and dipolar parts, respectively, of the Cr-Cr 
interaction. 

APR and EPR lineshapes have been measured for four transi- 
tions as a function of Cr concentration. APR second moments 
have been experimentally determined for the first time; APR 
and EPR half-widths have been compared. 

The EPR and APR half-widths quantitatively appear to be the 
same. 

At low concentration, the second moment is concentration 
independent and is much larger than the theoretical value 
predicted by Laurance, et pleas 

At high concentration, the APR second moment is much too 
small if all neighbor contributions are taken into account. 
We have shown that the nearer neighbor contributions to the 
second moment lie in the unmeasured wings of the line, thereby 
confirming previous predictions. A rough calculation has shown 
that qualitative agreement can be achieved when neighbors out 


to about 6A and 8A for the dipolar and exchange parts of the 


interactions, respectively, are excluded. 
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APPENDIX 


CRYSTAL BONDING 


Inasmuch as considerable time and effort were spent in perfecting 


a technique for bonding transducers, we present in this appendix our 


final method. We do not claim that this is the only successful approach, 


nor do we even claim that all of our steps are absolutely essential. 


We do, however, assert that the method was repeatedly successful. The 


procedure is set forth in step-by-step form: 


i 


Clean the quartz and ruby faces. The ruby could be moved with 
ordinary tweezers; however, in order to prevent breakage of 
the ultra-thin quartz plate, we made up a set of thin teflon 
tweezer tips for gentler handling. (A vacuum tweezer was 
also tried, but this was unsuccessful.) The cleaning solutions 
are: 
ay Concenbrawed nitric acid) min, 
b) Flush with distilled water. 
c) Detergent and water - scrub with Q-tip or ultrasonic 
cleaner for 5 min. If the quartz is placed on a wet 
tissue over a flat glass surface, it can be safely scrubbed 
with a Q-tip. An ultrasonic cleaner can break thin quartz 
plates. 
d) Fil@sh with distilled water. 
e) Reagent grade acetone - scrub for 5 min. 
Prepare the bonding material. We used Nonaq Stopcock Grease 
(Fisher Scientific Co.). Several bottles were stored in an 
ordinary refrigerator freezer at all times. Since Nonagq absorbs 


water from the air, the bottle should never be left open. 


or 


a) Place fresh blob of Nonaq on a glass slide and cover with 
watch glass. 

b) Put the slide on a small electric hotplate set to about 
210°C. Heat till granularity disappears - about 5 min. 
This step is absolutely essential. 

Just prior to bonding, give quartz and ruby one last wipe 

with reagent acetone. Make sure no drying marks remain on 

either crystal. 

With a glass rod, remove a small drop of Nonaq from the 

heated blob. Place drop on ruby face. Invert ruby and 

pick up quartz. 

Place quartz-ruby in a jig so that vertical pressure can be 

applied to the quartz. We used a staking tool; on the end 

of the plunger rod was a cylinder of silicon rubber, slightly 

smaller in diameter than the quartz. Press with about 150g 

for 15 min. Too long a time will cause vacancies in the bond 
material. The best bonds occur when the Nonaq spreads to edge 

Ofsaquartz but no further. 

Remove bond and check for uniformity by observing fringes, 

visible by an fluorescent desk lamp or, to a lesser extent, 

by incandescent light. 

Check on 500 MHz test rig. An even pulse envelope shape is 

indicative of good uniformity. The bond is now ready for use. 

Cooling down to low temperature must be slow in order to give 

the cold Nonaq a chance to flow. Do not evacuate the sample 

space, however: Nonaq tends to bubble. Our cooldown times 


were about 12 hours. 
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There is a large amount of art in successful bonding. We 
have had difficulty bonding on humid days, and sometimes a 
half-full bottle of Nonaq had to be discarded. On the other 
hand, what we thought were rather poor quality bonds occa- 
sionally gave excellent echoes. Trial and error sometimes 


seems to be the only answer. 
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